The excitation spectrum of mesoscopic proximity structures 
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We investigate one aspect of the proximity effect, viz., the local density of states of a superconductor- 
normal metal sandwich. In contrast to earlier work, we allow for the presence of an arbitrary 
concentration of impurities in the structure. The superconductor induces a gap in the normal 
metal spectrum that is proportional to the inverse of the elastic mean free path Zn for rather clean 
systems. For a mean free path much shorter than the thickness of the normal metal, we find a gap 
size proportional to Zn that approaches the behavior predicted by the Usadel equation (diffusive 
limit). We also discuss the influence of interface and surface roughness, the consequences of a non- 
ideal transmittivity of the interface, and the dependence of our results on the choice of the model 
of impurity scattering. 



I. INTRODUCTION 

A normal metal in good metallic contact to a supercon- 
ductor acquires superconducting properties like infinite 
conductance and the Meissner effect (seea and references 
therein). This so-called proximity effect has been exten- 
sively studied in the past decade mainly in the two limit- 
ing cases of fully ballistic propagation and diffusive mo- 
tion. Whereas the first may be realized experimentally 
in two-dimensional electron gases, the latter is realized 
in the most common samples made of semiconducting or 
structured metallic films. 

Although the ballistic and the diffusive cases provide 
useful bounds, real- world mesoscopic samples are often 
in the intermediate regime. For instance, the diamag- 
netic response of mesoscopic proximity cylinders (super- 
conducting wires covered by |-a|-,normal metal) h-as, at- 
tracted a lot of experimentalErcl and theoreticaltTO in- 
terest. It turned out that the experimental results could 
only be understood by considering intermediate impu- 
rity conceotrationaj. Previous-theoretical analyses in the 
clean limitH and the dirty limitM were not able to explain 
the experimental results quantitatively. Only recently it 
was shown that a qualitative different behavior emerges 
in the regime of intermediate impurity concentration^. 
With the help of this calculation a quantitative agrees 
ment with the experimental results could be obtainedcl 
and it was possible to determine the mean free path Zn- 

The precise determination of the degree of disorder is 
relevant for a characterization of the samples. It may 
serve as a basis for the understanding of the the exper- 
imental observation of a low-temperature paramagnetic 
reentrance effeciBS. This has already stimulated theo- 
retical suggestions that orbital currents might lead to a 
paramagnetic contribution0fi3. These currents depend 
on the degree of disorder in the normal metal and fur- 
ther investigations are necessary to quantify the influence 
of disorder. 

In the present work, we want to focus on another as- 
pect of the proximity effect, viz., on the change of the 
excitation spectrum of a normal metal connected to a 
superconductor. In particular, we are interested in the 
changes of the spectrum due to disorder. 



In this paper, we report on a comprehensive study on 
the density of states in a moderately dirty proximity sam- 
ple. The paper is organized as follows. In the next sec- 
tion we introduce the geometry and the parameters of 
our model. In Sec. Ill we present results for the density 
of states for different disorder concentrations and mod- 
els. Finally, we discuss the effect of rough surfaces and 
rough interfaces in Sec. IV. 



II. GEOMETRY AND MODEL 

The sample geometry that we have in mind is shown 
in Fig. |l|: we consider a slab geometry in which a normal 
metal layer of thickness d is connected by an ideal inter- 
face to a superconductor. The outer surface of the sand- 
wich is supposed to be specularly reflecting. In the fol- 
lowing, we will discuss the local density of states (LDOS) 
N(E, x) of this structure for a variety of physical situa- 
tions. The quasiclassical formulation of superconductiv- 
ity is most suitable to calculate the LDOS. To this end 
we have-soissed the real-time version of the Eilenberger 
equatiorjB'B 



- v F V.g(v F , r, E) = [-iEt 3 + A(r)fi 



a(r,E),g(v F ,r,E)} 
(1) 



for the quasiclassical 2x2 matrix Green's function g 
numerically (seeo and Appendix |A] for additional details 
of this method). Here and in the following, fi is set to 
one except for the final results. In Eq. (Q), A is the 
pair potential and a the impurity self-energy which has 
to be determined self-consistently by a scattering matrix 
equation. The Pauli matrices t% are used as a basis for 
the 2x2 matrix equation, vp is an arbitrarily oriented 
unit vector times the Fermi velocity. The quasiclassical 
Green's function is normalized according to 



g 2 (v F ,r,£0 = l 



(2) 



In principle, the pair potential A has to be determined 
self-consistently in the superconductor, whereas it van- 
ishes by definition in the normal metal. Since we are 
mainly interested in the properties of the normal metal 



1 



in the limit d 3> hvp/A, the spatial dependence of the 
pair potential in the superconductor plays no role. Thus, 
we approximate the pair potential by a step function 
A(r) = Ad(-x). 

The effect of impurities give rise to the self-energy er(r). 
In general the impurity self enetgy can be found from a 
solution of a t-matrix equationll3 In our particular case 
this equation can be solved and presented as 



a r 



E) 



n imp (r)v 2 (r)N (g(r, E)) 
l + v 2 (r)N 2 (g(r,E)) 2 



(3) 



This impurity self-energy contains two parameters, the 
impurity concentration ni mp and the strength of the scat- 
tering potential v(r) (that may be spatially dependent if 
there are different impurities in different parts of the sam- 
ple). Additionally, the normal metal-density of states at 
the Fermi level No enters into Eq. (|^). 

In the rest of the paper we will mostly use the Born 
approximation valid in the limit of weak scattering, where 
the denominator in Eq. (|J) can be neglected: 



a(v,E) 



2r imp (r) 



(4) 



Here, 



' imp 



we have introduced the elastic scattering time 
= (2n imp (r)w 2 (r)A r o) _1 which is related to the mean 



free path I — VF7i mp . Since we will also consider inhomo- 
geneous impurity distributions (on a mean level), we still 
allow for a spatial dependence of elastic scattering time 
and mean free path. 



ideally transmitting 
interface 



specularly reflecting 
surface 




non-magnetic impurities 



FIG. 1. Geometry of the SN interface. Zn (Is) is the elastic 
mean free path of the normal (superconducting) part. In the 
first part of the paper, the interface is supposed to be ideal 
and the surface specularly reflecting. 



Before we discuss our results, let us briefly summa- 
rize the known spectral properties in the limiting cases 
of clean and dirty limit. The case of a clean normal 
metal was already discussed in the 60's: the LDOS is in- 
dependent of the location in the normal metal, .vanishes 
at the Fermi energy and rises linearly close to itt§ll3. Its 
peculiar zig-zag from is shown in Fig. ^. The characteris- 
tic energy determining the jumps in the spectrum is the 
Andreev energy Ea- This characteristic energy follows 
from a semi-classical quantization condition for the An- 
dreev bound states. According to this rule we have to 
add the phases accumulated along a trajectory burying a 
bound electron-hole state. This phase difference is given 
by twice the phase shift of an Andreev reflection at the 
superconductor, which at E <C A is given by 7r/2. On 
the path through the normal metal both electron and 
hole accumulate an additional shift of 2 x 2Ed/vp, pro- 
portional to twice the time spent in the normal metal. 
Adding all contributions to the semiclassical phase and 
requiring it to be an integer multiple of 2tt leads to the 
characteristic energy of the lowest level 



E A 



Ad 



(6) 



The subsequent levels are (approximately) equidistant 
with level spacing 2Ea as can be seen from Fig. || (the 
deviations from the ideal level spacing are caused by the 
finite value of A * 5Ea)- 

^ / d = 0.25 



clean SN density of states (DOS) 
bulk normal metal DOS 
bulk superconductor DOS 




FIG. 2. Density of states of a clean normal metal connected 
to a superconductor. In the clean case, N(E,x) — N(E), 
i.e., there is no dependence on the spatial coordinate. The 
coupling to the superconductor leads to a linear suppression 
of the DOS at the Fermi energy. Here, we chose A = 5Ea- 



The set of equations (|l])-(0) has been solved numeri- 
cally. The integration oi-Eq. [l| was performed using the 
Riccati parametrizationtil (see Appendix |Al[ ) . The self- 
energies had to be determined self-consistently for all en- 
ergies in an iterative scheme. Finally, the local density 
of states is given by 

N(E, r) = N Rc(^Tvr 3 g(v, E)) . (5) 



Another well-known result on the spectrum has been 
obtained in the dirty (diffusive) limitliZHiS. In this case, 
the LDOS is characterized by a minigap in the spec- 
trum that is of the order of the Thouless energy Erh — 
TiD/d 2 , here D is the diffusion constant of the nor- 
mal metal and d its thickness. The LDOS of a meso- 
scopic superconductorpDormal metal sample was deter- 
mined experimentalrap in the presence of a magnetic 
field, and our theoryE3 led to a satisfactory understand- 
ing of those experimental results. 
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The induced minigap has also been discussed by field- 
theoretic means, seeE3. The relation of this gap to quan- 
tum chaos is discussed in RefsEi. 

A qualitative picture of the formation of the minigap 
follows from considering the Andreev bound states in the 
normal layer. An electron and hole traversing a diffusive 
trajectory in the time reversed direction are transformed 
into each other by Andreev reflection when hitting the 
superconductor. The coherent superposition of two sub- 
sequent reflections result in a bound state. The phase 
shift due to Andreev reflection is similar as in the clean 
limit discussed above. During the motion the two quasi- 
particle gain an additional phase shift 2ESt, where St 
is the time spent in the normal metal region and 2E is 
the energy difference of electron and hole. In a diffu- 
sive system the maximal time spent in the normal metal 
is ~ 2d 2 /D, i.e., twice the time to diffusion time. This 
upper time limit results in a lower bound to the energy 
~ ttD/AcP above which constructive interference is pos- 
sible. The minigap is thus expected to be approximately 
E g w irhD/Ad 2 = 0.785HD/d 2 . The numerically exact 
expression for the minigap in the dirty case is given by 
Eg = 0.780hD/d 2 , where the numerical factor is the so- 
lution of a transcendental equation (see Appendix |^). 
This result is in very good agreement with the estimate 
given before corroborating the simple picture of Andreev 
bound states. 



III. ARBITRARY IMPURITY CONCENTRATION 

How is the linear rise of the LDOS for the clean system 
transformed into the minigap in the diffusive system as a 
function of impurity concentration? To answer this ques- 
tion, we have solved the real-time Eilenberger equation 
including an impurity self-energy of the form a — (<?} /2r 
(Born approximation). The impurity self-energy was de- 
termined iiL | a_self-consistent way. 

We findEE3 that a gap forms at arbitrarily small im- 
purity concentrations. This is shown in Fig. |^: even for 
values of the elastic mean free path Zn that are 30 times 
larger than the normal-layer thickness, the formation of 
the low-energy gap is clearly visible. The gap increases 
with 1/Znj saturates for Zn ~ d and then decreases again 
as expected from the dirty- limit theory since D ~ vfZn, 
see Fig. ^. The gap does not depend on the location in 
the normal metal as can be seen in Fig. i.e., it is a 
global feature. The shape of the LDOS, i.e., its depen- 
dence on energy, however, varies on traversing the normal 
layer. 

The existence of a minigap is in-Jine with qualitative 
considerations given by McMillanEl He argued that - 
quite generally - the density of states should show a 
gap of order of the inverse of the escape time, i.e., the 
time an electron spends in the normal layer before being 
Andreev-reflected. If we replace the escape time by the 
scattering time in the (almost) clean system and by the 



diffusive escape time in the dirty system, we obtain the 
non-monotonous behavior shown in Fig. y. 



- x/d = 1 .0 surface 

- x/d = 0.8 
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FIG. 3. LDOS N(E, x) for different ratios of mean free 
path Zn to normal-layer thickness d. The minigap is constant 
throughout the normal metal, but the energy dependence of 
the LDOS changes with location. 



i; /d = 0.01 




FIG. 4. Size of the minigap obtained by numerically solving 
the Eilenberger equation. 



We have also replaced the Born model for impurity 
scattering by the unitary limit which follows from the 
t-matrix approximation for s-wave scattering Eq. (^) in 
the limit vNq 3> 1. To facilitate a comparison we used 
for the parameter m/No, characterizing the strength of 
the impurity self-energy in the unitary limit, the value of 
l/2r from the Born approximation. Figure [| shows the 
minigap for both the Born and the unitary approxima- 
tion. The minigap is slightly reduced in the unitary limit, 
but its functional dependence on the mean free path is 
practically unchanged. 
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FIG. 5. Dependence of the induced minigap on bulk dis- 
order for two models of impurity scattering, viz., Born and 
t-matrix approximation. The differences are insignificant, i.e., 
the minigap is stable and is not dependent on the choice of 
the Born approximation. 



IV. INTERFACE AND SURFACE ROUGHNESS 

Real-life interfaces and surfaces of typical proximity 
samples are rough (for afi example, see the photograph in 
the second paper of RefE). In the quasiclassical language 
such a roughness will lead to a mixture of different tra- 
jectories, thus smearing out the singular Andreev bound 
states on a given path. A convenient way to include 
this effect in the quasiclassical formalism is to model the 
rough surface by a thin dirty layer which covers the in- 
ner side of the surface, see Fig. g. The thickness 8 of the 
layer should be taken to be much smaller than all char- 
acteristic lengths of the rest of the system (in in our case 
this are the thickness of the normal metal d and its mean 
free path The disorder in the layer is included by 
an impurity self-energy with a mean free path ii aye r- At 
the outer surface inside the dirty layer specular reflec- 
tion is assumed and the Green's functions in the dirty 
layer are continuously connected to the ones in the nor- 
mal metal. Under these conditions the ratio l\ aye -[/8 is the 
only parameter measuring the roughness of the real-life 
interface. 

Similarly the interface roughness between the normal 
metal and the superconductor is modeled by a thin dirty 
layer residing now at the interface. The layer is charac- 
terized by the same parameters as before. Now we have 
continuity of the Green's function on both sides of the 
dirty layer. 

It is clear that this model for roughness should be taken 
with care, since it is by no means microscopically justi- 
fied. The parameters characterizing the dirty layer are 
not related to the real parameters of the interface. It is 
however clear that this model covers the essential prop- 
erties of a realistic rough surface, i.e. it couples classical 
trajectories which would be uncoupled for a specular in- 
terface. Thus it leads to an averaging over spectral quan- 
tities from different regions of the system in the same 
way as a realistic system does on average. One should 
however keep in mind, that we can only determine aver- 
ages over many characteristic lengths on the roughness 



in this way. E. g. in small cavities with a characteris- 
tic length of roughness which is comparable to the size 
of the system, such an averaging is not appropriate and 
fluctuations may become important. 
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FIG. 6. We model a rough surface (or interface) by replac- 
ing it by a thin dirty layer. The scattering time r in this layer 
and its thickness are chosen such as to lead to the proper scat- 
tering characteristics of a rough surface. 



A. Surface disorder 

The result of a calculation for a rough outer surface 
is shown in the left plot of Fig. [j]. We have kept the 
roughness of the surface fixed in this plot and varied 
the bulk disorder. Only the LDOS at the outer surface 
is plotted. In our example, the roughness parameter is 
such that there are on average two scattering events in 
the dirty layer, i.e., the outer surface is definitely non- 
specular. The induced minigap is not very susceptible to 
the presence of the surface roughness as follows from a 
comparison with Fig. |^. Using McMillan's argument that 
the minigap should be inversely proportional to the es- 
cape time, we can qualitatively understand this behavior. 
The effect of the scattering by the rough surface leads to 
a reduction of the number of the shorter trajectories, but 
the length of the shortest trajectory itself is not changed. 
This can be most clearly seen for an intermediate bulk 
disorder / w 3d, where the sharp increase at the mini- 
gap for the specular surface is nearly absent. The rough 
surface also leads to a change in the spatial form of the 
density of states, which is not shown in the plot. 

We remark in passing that surface roughness without 
bulk disorder will not lead to the formation of a mini- 
gap, since there will be no upper cutoff for the trajectory 
lengths in this case. Nevertheless it influences the energy 
dependence of the LDOS, by virtue of a similar effect as 
mentioned previously. Due to the reduction of the num- 
ber of the longer trajectories spectral weight is shifted 
to higher energies. This can already be seen from the 
solid curve in the left plot Fig. 0, where the linear energy 
dependence in the case of specular reflection has turned 
into some weaker energy dependence. 
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B. Interface disorder 



decreased. 



In many experimental situation the interface between 
a normal metal and a superconductor will be nonideal in 
the sense that either transport through the interface is 
not along classical trajectories or the electrons are par- 
tially reflected. This can be due to either an oxide layer 
or an alloy in the interface region, and/or due to dif- 
ferences in effective masses or Fermi velocities between 
the superconductor and the normal metal. The generic 
properties of these disordered interfaces can be included 
in the quasiclassical formalism by a thin dirty layer lo- 
cated at the interface. This model will also include the 
effect of a finite reflection at the interface. Since it is a 
well know universal property of disordered conductors to 
have a bimodal distribution of transmission eigenvalues 
(see e. g. 24), we believe that the model of such an disor- 
dered layer will cover most of the characteristics of rough 
contacts. 




FIG. 7. LDOS at the outer surface (a; = d) under the 
influence of surface disorder (left plot) or interface disorder 
between superconductor and normal metal (right plot). The 
thickness of the layer is chosen to be 8 = 10 _5 d and the mean 
free path in the layer is Zi aycr = 8/2. The different curves in 
each plot are for different degrees of disorder in the normal 
metal as indicated in the right plot. A comparison with Fig. |^ 
shows that the size of the minigap is practically conserved 
in the left plot, whereas the non-ideal interface leads to a 
decreased minigap in the right plot. 



In the right plot of Fig. we show the effect of a thin 
dirty layer at the interface between superconductor and 
normal metal. The interface is chosen to have a mean free 
path Zn = 5/2, similar to the previous case of the rough 
surface. This is supposed to mimic a rough interface 
between the metals with equal Fermi velocities. We find 
a significant reduction of the minigap. This is clearly seen 
from comparison with the left plot of Fig. |t] and shown as 
a function of the interface roughness parameter in Fig. ^. 
For roughness parameters 5/l\ ayeT smaller than w 1 the 
reduction of the minigap is weak. If this parameter is 
of the order of or larger than 1 the minigap is strongly 
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FIG. 8. Influence of interface roughness (parametrized by 
8/ha.yev) on the induced minigap in the normal layer. In con- 
trast to surface roughness, interface roughness leads to a pro- 
nounced suppression of the minigap. In particular roughness 
parameters larger than 1 lead to a strong suppression of the 
minigap. 



This behavior can be understood by looking at the 
reflection properties of the rough interface. If the rough- 
ness parameter is smaller than 1 the reflectivity will be 
roughly given by R « 5/li aycr , i.e., the probability for a 
particle to be reflected is given by the ratio of the mean 
time spent in the layer to the scattering time. As a rough 
estimate of the effect of the finite reflectivity R > for 
the minigap, we use the fact that the mean time that an 
electron spends in the normal metal will be increased by 
a factor 1 + R+R 2 + .. = 1/T, where T is the transmission 
eigenvalue of the interface. Hence, the minigap will be 
reduced by a factor T w (1 — 5/l\ aycl ). The influence of 
finite reflection coefficients in the clean system was also 
discussed inE£l. 

If the interface is strongly disordered, i.e., the rough- 
ness parameter is larger than 1, it behaves more like a 
disordered metal, having a bimodal distribution of trans- 
mission eigenvalues. Most of the transmission eigenval- 
ues are close to zero and a few are close to unity al- 
lowing for Andreev reflection. The ratio of the number 
of closed to the number of open trajectories is roughly 
given by Thus, an electron typically is normally 

reflected by the dirty layer many times until it has the 
chance to be Andreev reflected. This strongly enhances 
the length of the trajectories and therefore reduces the 
minigap roughly by a factor Zi ayer /<5. This explains the 
strong suppression of the minigap that sets in for values 
of the roughness parameter larger than 1 as shown in 
Fig. |. 

A reduction of the minigap appears in the same man- 
ner, if the superconductor is dirty, i.e., has an elastic 
mean free path lg shorter than the coherence length 
Co = wf/As- The disorder close to the interface acts 
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similarly as the dirty layer at the interface. The role of 
the "roughness" -parameter is now taken over by £oAs> 
since Andreev reflection occurs in this layer of the super- 
conductor. Again the effective length of the trajectory 
in the normal layer is enhanced by normal reflection at 
the impurities in the superconductor. Thus, we expect 
a qualitative similar behavior that the minigap will be 
reduced by an increase of the disorder in the supercon- 
ductor. 



V. CONCLUDING REMARKS 

In conclusion, we have studied the local density of 
states of a proximity sandwich for a variety of situations. 
At arbitrarily small impurity concentration, a gap opens 
at the Fermi energy; it is maximal if the elastic mean free 
path is of the order of the normal-layer thickness. We 
have numerically calculated this gap and its dependence 
on surface and interface roughness. Whereas the gap is 
relatively stable to surface roughness, it is strongly sup- 
pressed by interface roughness. We have also investigated 
the effect of a non-ideally transmitting interface. Lastly, 
we have investigated the influence of different models of 
impurity scattering (Born vs. unitary limit) and shown 
that the two models lead to a qualitatively similar be- 
havior. 
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APPENDIX A: NUMERICAL CALCULATIONS 

For our calculations we use the Riccati parametrization 
of the Eilenberger equations (seeE-3). We represent the 
Green's function on a trajectory in the form 



1 
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1 — aa^ 2a 
aat V 2cjt aa^ - 1 



(Al) 



We have introduced the Andreev amplitudes a and (V, 
which depend on the variables as g. From the Eilenberger 
equation ([!]) one derives two decoupled numerically sta- 
ble equations of the Riccati type: 



-vpVo = 2Cja -\ 
vpVat = 2uia' f 



A* a 2 - A , 
f Aa t2 - A* 



(A2) 



Here, the impurity self-energy is included in u> and A: 
ui = —iE + (Tn(r) 



A = A(r) + ffi 2 (r) 



(A3) 



We note that the first equation is stable in an integra- 
tion in positive direction along the trajectory, whereas 
the second equation is stable for integration in the op- 
posite dir ecti on. In the stable direction the differential 
equation (A2) is conveniently integrated by a discretiza- 
tion which leads to the very accurate expression 



A - (f2 + Q)a n ) e - 2J W«F - A-(n- Cu)a 7 , 



(A4) 



Here, fl = (A 2 + a) 2 ) 1 / 2 , and h is the step size. 



APPENDIX B: MINIGAP IN THE DIRTY LIMIT 

The Usadel equation! 2 ! 

D d 2 



2 dx 2 



e 



' sin — A cos 



(Bl) 



contains only one energy scale, the Thouless energy 
Et\i — fiD/d 2 (apart from the pair potential A). The 
minigap that emerges on this scale can be explained 
by the following argument. The LDOS is given by the 
real part of g = cosG. If it vanishes, we may write 
cos 8 = — isinht? where i) is real, and the Usadel equa- 
tion in the normal metal reduces to 



de 



2Ed 2 
D 



cosh d 



(B2) 



(here, £ = x/d). This differential equation is of elliptic 
type and can be integrated twice. Next, we apply the 
boundary conditions $(£ = 0) = and ^^(£ = I ) = 0. 
The latter accounts for the symmetry of the reflecting 
surface. The first condition approximates $ by its value 
in the bulk superconductor at low energies which is jus- 
tified, if A is sufficiently large. We find that Eq. ( |B2] ) is 
only solvable, if 



E > 0.780hD/d 2 



(B3) 



The minigap of the dirty limit is set by this restriction. 
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